Critical dimensions of the diffusion equation 
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We study the evolution of a random initial field under pure diffusion in various space dimensions. 
From numerical calculations we find that the persistence properties of the system show sharp tran- 
sitions at critical dimensions d\ ~ 26 and d2 — 46. We also give refined measurements of the 
persistence exponents for low dimensions. 
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Diffusion is a ubiquitous process in Nature and is of 
fundamental importance in physical, chemical, and bi- 
ological systems. Pure diffusion, as described by the 
simple heat equation, has been so well-studied over so 
many years that it is difficult to imagine there is still 
much about it to be understood. Recently, the persis- 
tence probability, g(i), of the diffusion equation in one, 
two, and three dimensions was measured and found to 
be described by non-trivial power laws ||l|]: q{t) ~ i~^('^\ 
where d is the dimension of space. As yet there is no 
exact prediction for the exponents 9{d). 

The study of persistence is a rapidly emerging field [^ , 
the purpose of which is to probe the statistical nature 
of dynamical processes with more resolution than that 
allowed by conventional measures, such as autocorrela- 
tion functions, ft has been applied to a wide range of 
systems including phase ordering, interface kinetics, and 
the diffusion equation. Direct experimental observation 
of persistence has been achieved in the phase ordering 
of liquid crystals ^ and just recently by NMR studies 
of the diffusing polarization of a one dimensional gas [Q . 
Strictly speaking, the persistence probability measures 
the probability that a particular event has not occurred 
up to time t. More loosely, by "persistent properties" , 
we include other statistical measures, such as the sign- 
time distribution (STD) [|,§. The STD is a histogram 
measuring the fraction of time t a bimodal variable has 
been in one of its two states. Examples are the fraction 
of time an Ising spin is "up" , or the fraction of time a 
stochastic variable is above its mean. The tails of the 
STD implicitly encode the persistence probability. The 
shape of the STD yields much useful information about 
the spatio-temporal structure of the dynamical process. 
For example, it can be used to infer long-lived features 
in interfaces undergoing kinetic roughening M. 

In this Letter we introduce a simple algorithm which 
allows us to study the persistence properties of the diffu- 
sion equation in arbitrary dimensions. Our main result is 
that there exist two critical dimensions in the behavior of 
the STD. We find that for dimensions below di ~ 26 the 
STD is purely concave (and thus typical sign-time histo- 
ries dominate q{t)), while for dimensions above ^2 ~ 46 
the STD is purely convex (and thus q{t) is controlled by 
rare sign-time histories). For dimensions in the range 



di < d < d2 the STD has more than one extremum and 
the description of q{t) in terms of typical and rare events 
is less straightforward. 

We consider a scalar field (?!)(x, t) in a d-dimensional 
space which is to be evolved under the diffusion equation: 



dt 



DV" 



(1) 



where we have (/>(x, i = 0) = V'(^)- The random field ip is 
taken to have zero mean and to be Gaussian distributed 
with correlations {ip{x)ip{x')) = A(5''(x — x'). The results 
that follow are completely independent of the nature of 
the distribution (so long as it describes short-range cor- 
relations) and the values of D and A (so long as they 
are positive definite). We define the persistence proba- 
bility for this system as the probability that the diffusion 
field at a given point in space (the origin, say) has never 
changed sign up to time t. The sign-time T(i) is simi- 
larly defined as the proportion of time t the field at the 
origin has been positive. The STD is the histogram con- 
structed from the ensemble of sign-times collected from 
a large number of realizations. 

The persistence probability q{t) for the diffusion equa- 
tion has been measured previously from explicit sim- 
ulations of a discrete diffusion process on a finite d- 
dimensional grid pi for dimensions d=l, 2, and 3. For 
a given system one collects the persistence time for each 
site of the lattice and forms q{t) by binning these values. 
There is an implicit assumption concerning self-averaging 
since these persistence times are not independent. This 
is equivalent to stating that q{t) also measures the frac- 
tion of persistent sites in the system, which is certainly 
plausible in the large time limit. 

We have measured q{t) and the STD using an alter- 
native algorithm which has many advantages over direct 
simulation of the discrete diffusion process. The algo- 
rithm is constructed as follows. First we note that the 
diffusion equation may be explicitly integrated. Restrict- 
ing our attention to the diffusion field at the origin we 
have 



(2) 



0,t)= / d^y giy,t)iP{y) 



where g(x,t) = {4iTDt)~'^/'^ exp{-x'^ /ADt) is the heat 
kernel. The key observation is that the solid angle in- 



tegration in Eq.(g) may be absorbed into a redefinition 
of the random field, as originally noted by Hilhorst [Q . 
Defining a Gaussian random field ^(r) with correlations 

{^[r)'^[r')) = d{r - r'), we have 



(0,t) (X dr r 



Jd-l)/2-r 



'/**M , 



(3) 



where we have scaled out D and A and ignored the (pos- 
itive definite) prefactor, as only the sign of (l>{0,t) is rel- 
evant to the persistent properties. To measure q(t) for a 
given dimension d we construct a random field \l/(r) and 
perform the radial integral above for a series of incre- 
mented times, until the integral changes sign. We record 
the time at which this sign change occurs and repeat the 
procedure for another realization of ^(r). This collection 
of persistence times is then used to construct q{t). If we 
wish to measure the STD we evaluate the integral for a 
series of incremented times within a fixed time range and 
record the fraction of time for which the integral is posi- 
tive. The computer time required for these procedures is 
essentially independent of the value of d. 

This algorithm is superior to direct numerical simula- 
tion in many respects: i) any dimension may be studied 
(even non- integer), ii) there are no finite size effects, iii) 
there is no implicit assumption concerning self-averaging, 
iv) in measuring q{t) a given realization is discarded as 
soon as the field changes sign (in a direct simulation, all 
field points are studied even if the vast majority are no 
longer persistent), v) the algorithm may be optimized 
allowing very long time intervals to be studied with ease. 

We shall now discuss a few technical details concerning 
optimization. One may think that in order to measure 
q{t) one must evaluate the integral in (0) on a fine grid of 
incremented times, so as not to miss any event in which 
the field changes sign twice over a short time interval. 
In fact this is not the case. For a given realization of 
^ we have explicitly evaluated the function (j>{0,t) on a 
fine mesh of time slices and find it is smooth, becoming 
more so as time proceeds. This is reasonable when one 
considers that under diffusion the field at a given point 
will have increasingly smoother dynamics with time. Due 
to the progressive smoothness of (/)(0,i) we need only 
sample it on a logarithmic time scale. 

The base of the logarithm to be used is determined 
from calibration. We find the exact zeroes of 10^ realiza- 
tions over three decades of time using a fine linear mesh 
of time slices. We then run the optimized code which 
searches for these zeroes first using a logarithmic scale, 
and then using recursive bisection of a logarithmic inter- 
val if it has reason to believe a zero lies within. The loga- 
rithm base is chosen so that all of the zeroes from the 10^ 
samples are found correctly. This limits the systematic 
error of our results to less than 0.01%. The core of the 
algorithm concerns the detection of zeroes within a loga- 
rithmic interval, based on the value of the function 0(0, t) 



at the boundaries of the interval. We use two techniques. 
The faster, but less precise, algorithm uses recursive bi- 
section of an interval only if the values of 0(0, t) at the 
boundaries have different signs (in which case there has 
to be an odd number of zeroes in the interval). The more 
sophisticated algorithm also uses information about the 
derivatives of 0(0, t) at the boundaries of the interval. 
This algorithm will also use recursive bisection to check 
for a zero if the values of 0(0, t) at the boundaries have 
the same sign, but the derivatives have opposite signs 
(indicating an extremum in the interval and the possi- 
bility of an even number of zeroes). As an example of 
calibration, for d = 1 the first algorithm finds all zeroes 
correctly over three time decades if the logarithm base 
used is 1.1. The second algorithm finds all zeroes cor- 
rectly with the logarithm base 2.0. (As d increases the 
second algorithm requires a smaller logarithmic base to 
find all the zeroes - as an example, for d = 20 the re- 
quired base is 1.2.) For measurements of q{t) we always 
use the second algorithm, but for measuring the STD we 
generally use the first algorithm (as the STD requires 
much longer computing time). 

Using this method we have been able to measure q(t) 
and the STD for a wide range of dimensions. For dimen- 
sions d=l, 2, and 3 we are able to make estimates for the 
persistence exponent 9 as shown in Table 1, with much 
higher precision than formerly obtained H. In Table 1 
we also give estimates for 6 in higher dimensions, which 
until now have not been measured. By comparing these 
results with those obtained from the independent interval 
approximation (IIA) m we see that the IIA consistently 
underestimates the exponent values. 



d 


m 


d 


9{d) 


1.0 


0.12050(5) 


10.0 


0.4587(2) 


2.0 


0.1875(1) 


20.0 


0.6556(2) 


3.0 


0.2382(1) 


30.0 


0.8053(3) 


4.0 


0.2806(2) 


40.0 


0.9312(4) 


5.0 


0.3173(2) 


50.0 


1.0415(5) 


26.0 


0.7491(3) 


46.0 


1.0010(5) 



Table 1: Numerically calculated values of the persistence 
exponent for a range of dimensions. 

In Figures 1 and 2 we present log- log plots of q{t) for a 
sample of dimensions to show the quality of our data. We 
generally calculate 0(0, t) over six decades of time, and 
average over 10^ realizations. Our exponent estimates 
were made in the same fashion as in Ref. M, namely, by 
plotting t^q{t) versus \ogt and adjusting the exponent 
until the data is as level as possible. The error bars are 
estimated from the resolution of this procedure. 

For the diffusion equation it is known that for large d 
the persistence exponent diverges as ^/d m. By plotting 
our measured exponents on a log-log scale we have fitted 
this behavior and find 6 -- c\/d with c = 0.1475(10). At 



the bottom of Table 1 we have given exponent estimates 
for the two critical dimensions in this system. To clarify 
the role of these dimensions we shall now discuss our data 
for the STD's. 




10 l(f- l(fi 104 10^ lO*' t 10'^ 
1. Persistence probability q{t) versus t, for dimensions 



FIG 
(from top to bottom) 1 



2, 3, 4, and 5. 




10 l(f- 10^ 10^ 10^ 10^ t 10'^ 
FIG. 2. Persistence probability q{t) versus t, for dimensions 
(from top to bottom) 10, 20, 30, 40, and 50. 

To measure the STD, P(r, t), we evolve the system for 
three decades of time and by locating zeroes of (f>{0,t) 
determine the sign-time r(t) for a particular realization. 
This process is repeated for 10^ samples and the STD 
is constructed by binning the sign-times. In Ref. Iq] 
it is shown that the STD has the exact scaling form 
P{T,t) = (l/i)/(u) where u — r/t. In what follows we 
shall describe the STD in terms of the scaling function 
f{u), which is symmetric about u = 1/2. The tails of 
the STD encode information about the persistence prob- 
ability 0. For example, the left tail is expected to vary 
as u^~^. Thus, ioi 9 < 1 the STD will have integrably 
divergent tails, whereas ioi 9 > 1 the STD will have van- 
ishing tails. Since 9 diverges with d, there exists a critical 
dimension d2 at which 9{d2) = 1, and which separates 
(crudely speaking) convex and concave STD's. From our 
numerical work we find that this scenario is too simple, 
and that for a range of dimensions the STD has more 
than one extremum (as shown in Figure 3). 
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FIG. 3. STD's for dimensions 10, 20, 30, 40, and 50. 

We start by numerically determining d2- By simply 
measuring the persistence probability q{t) for a range of 
dimensions we can identify that dimension at which 9 = 
1. As can be seen from Table 1, we find that 6'(46.0) = 
1.0010(5). A closer analysis leads us to the result ^2 — 
45.9(1). This may be verified (to a lesser precision) by 
studying the STD's. In Figure 4 we show the tails of the 
STD's for d = 45.0 and d = 46.0. Note in the former case 
the upturn signalling a (very) weak integrable singularity 
in the distribution. 
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FIG. 4. Tails of the STD for d = 45.0 and d 
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: 46.0. 



From the STD's in Figure 3 we see there is another 
critical dimension di separating STD's which are purely 
concave from those with three extrema. There is no ob- 
vious way to identify this by measuring the persistence 
exponent. Instead we must use direct measurements of 
the STD's which naturally leads to a less precise esti- 
mate. In Figure 5 we show the measured STD's for a 
small range of dimensions around d = 26.0. From this 
data we arrive at the result di = 26.0(5). In an earlier 
work |q| d2 was estimated by assuming that the STD in 
d2 dimensions is perfectly fiat. This led to the predic- 
tion ^2 — 35.967 . . ., which, within the error bars, is the 
arithmetic mean of the measured values of di and ^2 ■ 

Our main result with regard to the persistence prop- 
erties of the diffusion equation concerns the non-trivial 
behavior of the STD as a function of d. For d < di = 
26.0(5) the STD is a purely concave function, and thus 
the persistence probability (which is determined by the 



tails of the STD) is dominated by typical sign-time his- 
tories. For d> d2 = 45.9(1) the STD is a purely convex 
function, and in this case the persistence probability is 
determined from rare events in the set of sign-time his- 
tories. There is a range of dimensions di < d < 6,2 for 
which the STD has three extrema (as shown in Figure 
3). In this case the persistence probability is still de- 
scribed by the most typical sign-time histories, but we 
see that sign-time histories around u = 1/2 (meaning the 
field (/)(0, t) spends half of its time above the mean) are 
relatively typical also, with the least typical events oc- 
curring for some d-dependent value of u. These results 
are universal. The STD is independent of details of the 
initial distribution of the diffusion field, and the value of 
the diffusion constant. It is also asymptotically (i.e. for 
large times) independent of higher-order processes, such 
as the corrections due to an underlying lattice. 
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FIG. 5. STD's for a small range of dimensions about 
d\ ~ 26.0. The horizontal lines are a guide to the eye. 

Clearly if one is to find applications of our high- 
dimensional results, then the space in which the diffusion 
process occurs cannot be taken as the familiar physical 
space (at least within the realm of classical physics). A 
very common manifestation of the diffusion equation is 
within the field of stochastic processes. Indeed, the one 
dimensional diffusion equation is the simplest possible 
form of the Fokker- Planck equation, describing the prob- 
ability distribution of a single random walker |10|. For 
N independent random walkers in one dimension, the 
multi-variate probability distribution p{xi, . . . ,xm) sat- 
isfies the A^-dimensional diffusion equation. Given ran- 
dom initial positions of the walkers, we should focus on 
the sign-time histories of p with respect to its mean value 
over the ensemble of initial conditions. The distribution 
of these sign times will have sharp changes in behavior 
as the number of walkers N is varied through A'^i = 26 
and N2 = 46. If each walker exists in a d dimensional 
space, then the effective dimensionality of the diffusion 
equation is Nd and the critical walker numbers A^i and 
N2 will be changed accordingly. Analogous statements 
may be made for a system of N directed polymers (e.g. 



magnetic flux lines) in a d dimensional space, since the 
multi-variate partition function for N lines also satisfies 
the Nd dimensional diffusion equation ||l^ . 

Our results also indicate that persistent fluctuations 
in the phase space of few-body systems may have sharp 
transitions on varying the number of degrees of free- 
dom. By "few-body" systems, we have in mind organic 
molecules composed of tens of atoms, the dynamical 
properties of which are currently receiving a great deal 
of attention |13] . 

In this Letter we have presented results for the persis- 
tence probability and sign-time distribution for the dif- 
fusion equation over a wide range of dimensions. These 
results have been obtained via a new algorithm which al- 
lows extremely precise measurements of the persistence 
properties. This algorithm is based on the integrability of 
the diffusion equation, a feature shared by many other in- 
teresting problems such as the Edwards- Wilkinson model 
of interface growth ||l^ , and the Burgers equation jlj] of 
fluid turbulence. The persistence properties of these and 
similar models may be analyzed to high precision using 
the methods outlined here. 
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